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Abstract 

The recently observed properties of quasiparticle resonant states near impu- 
rities on the surface of superconducting Bi2Sr2CaCu208+5 demonstrate tliat 
in-plane Andreev processes are either absent or phase-incoherent. Analysis 
of the spectral and spatial details of the electronic structure near a Zn impu- 
rity also suggest an effective magnetic component of the impurity potential. 
Further experiments are proposed to clarify whether the effective moments of 
nearby impurities are correlated. 



Over the past few years several authors have emphasized the wealth of information avail- 
able from local probes of impurity properties in correlated electron systems, and particularly 
in superconductors whose homogeneous order parameters (OP) are anisotropic in momen- 
tum [|l|. A parallel improvement in scanning tunneling spectroscopy (STS) has allowed this 
vision to become a reality through direct observation of the local density of states (LDOS), 
first in niobium 0, which has a momentum- independent OP, and this year in the high- 
temperature superconductor Bi2Sr2CaCu208+5 (BSCCO) which has an anisotropic 
OP. The electronic structure of the BSCCO surface is much more complex than that of the 
niobium surface; there are local moments in the copper-oxygen planes and, under certain 
conditions, a pseudogap state. Recent work ^ has emphasized the role of the pseudogap 
state in determining the properties of clean surfaces of high-temperature superconduc- 
tors at temperatures near Tc, and STS and photoemission have directly demonstrated its 
existence on the surface of BSCCO. The pseudogap state is characterized by a single-particle 
gap of dx2_y2 symmetry, but the Andreev processes one would expect in a superconducting 
state are (according to differences in mechanism) either absent or phase-incoherent. 

Here the STS measurements near impurities will be shown to unambiguously 

demonstrate essentially complete suppression of Andreev processes in the vicinity of the 
impurity, even at very low temperatures. This low-temperature indication of a pseudogap 
implies its importance to the nature of the superconducting state and the operation of de- 
vices (such as Josephson junctions) with such materials. Indications of an effective magnetic 
component can also be seen in Ref . . Thus STS provides a direct probe of both the lo- 
cal Andreev (superconducting) coherence and the local magnetic properties on the BSCCO 
surface. 

A brief review of the experimental results from Refs. is in order. Theoretical 

predictions which have been confirmed include the presence of quasiparticle resonances near 
nonmagnetic impurities in anisotropic OP superconductors [jlO|, as well as the suppression of 
the gap feature near the impurity and the asymmetry of the resonance peak due to the energy 



dependence of the quasiparticle density of states |Tl[]. The disagreements with previous 
theory, however, are striking. The most noticeable one is that the resonant state has only 
been detected on the hole side of the spectrum, both on the impurity site and everywhere 
else around the impurity. Whereas previous theories are consistent with an LDOS measured 
at the impurity which is entirely hole-like, these same theories unambiguously predict the 
LDOS at nearest-neighbor states will be almost entirely electron-like. Indeed these theories 
predict the spatially-integrated LDOS (or DOS) will be nearly particle-hole symmetric even 
though the LDOS at any particular site is not. A second unexpected element in the data is 
the presence of a second, much smaller, spectral peak on the hole side in Ref. [Q. Whereas 
previous theories are consistent with two resonances at a single impurity, the spatially- 
integrated amplitude of each resonance should be approximately the same, unlike what is 
seen in Ref. 0. 

The above two issues are disagreements between the experimental and theoretical DOS, 
but there are also two significant disagreements in the LDOS. The resonance has a large 
amplitude at the impurity site, whereas calculations indicate that the largest amplitude in 
the LDOS should occur at the nearest-neighbor sites. Finally, the gap feature is seen on the 
impurity site, where it does not appear in calculations. |l| 

The calculations here incorporate the pseudogap state into the evaluation of the dif- 
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ferential conductance (dl/dV). After evaluating the LDOS for several impurity potential 
models, allowing for spatial extent and magnetic character in the potential, the impurities 
of Refs 1^-0 are found to be highly localized and the four discrepancies above can be rec- 
onciled. Finally, in the case of Ref. a magnetic component to the effective Zn impurity 



potential is evident. [12 



The calculations of the LDOS are based on the Hamiltonian 

which includes a site-dependent potential which can be magnetic (V5), nonmagnetic (Vq), 
or a combination of both, i and j label sites and a labels spin. The homogeneous electronic 
structure, expressed as hopping matrix elements {tij) for the first five nearest neighbors, is 



taken from a single-band parametrization of photoemission data. ||T3[ Large variations in 
hopping matrix elements (> 50 meV) produce results very much at odds with experiment, 
whereas smaller variations are mimicked by the site-dependent potential. Hence such changes 
will be ignored here. Only on-site and nearest-neighbor order parameters A^- are nonzero, 
and the maximum OP on the Fermi surface, A^ax = 40 meV. 



The electronic structure of the inhomogeneous system (including the impurity) is deter- 
mined by direct numerical solution in real space of the Gor'kov equation (in Nambu form) 
for the inhomogeneous Green's function, G = g + gVG = (I — gV)~^g, within a real-space 
region around the impurity beyond which the potential is negligible. The Aj^'s are found 
self-consistently in this process, for they determine the off-diagonal components of the po- 
tential V. Spectra outside this real-space region are constructed according to the generalized 
T- matrix equation: G = g + gV [I — GV]~^ g. Once G has been calculated throughout the 
region near the impurity, the LDOS and DOS are obtained from its imaginary part and the 
lattice Wannier functions. Then 
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where 0tr(x; i) is the overlap of the Wannier function at site i and spin a with the STM tip 
at x, and nsTui!^) is the occupation function of the STM tip. Resonances correspond to 
new peaks in the differential DOS (the difference between the inhomogeneous and homoge- 
neous DOS); their energies are shown in Fig. [^(ab) for magnetic and nonmagnetic single-site 
impurities. 

If Andreev processes are suppressed, either by reduction in their amplitude or phase 
coherence, a resonance's DOS will become more electron-like or more hole-like. Reduction 
of the amplitude of the homogeneous anomalous Green's function f{i,j;uj), due perhaps 
to a local antiferromagnetic (AF) order, decreases the mean-field coupling between electron 
and hole excitations. Note that this is very different from the fully electron or hole-like 
character of the LDOS at the impurity, which originates from a vanishing f{i,i;uj) in the 
dx2-y2 state. Figure |l|(cd) shows the DOS of a resonance for three systems with a 40 meV 
dx2-y2 gap: a fully superconducting gap (solid line), a gap with a 25 meV superconducting 
component (dashed line), and a pseudogap with no superconducting component (dot-dashed 
line). As the superconducting component is reduced, the electron-hole symmetry diminishes. 
The nonmagnetic potentials of Fig. |l|(c) are chosen (1.375 eV, 1.000 eV, and 0.833 eV, 



3 



respectively) so the resonance peak is at —1.5 mV (the same as Ref. 0). The magnetic 
potentials in Fig. ||(d) are the same as those in Fig. ||(c). 

The reduction of the electron-like peak from phase decoherence is similar to the effect of 
amplitude suppression of f{i, i; u). For a resonance at u, the peak at —u comes from terms 
with products of pairs of anomalous Green's functions. For a phase incoherent pseudogap 
state the expectation value of these pairs, and thus the amplitude of the electron-like peak, 
is diminished. The effect of this is shown in Fig. ^(ef ) as a dashed line corresponding to 
partial (half) and a dot-dashed line corresponding to no phase coherence. The nonmagnetic 
potential is 1.375 eV in Fig. |I](e) and the magnetic potential is 1.375 eV in Fig. ^(f). Note 
that for the purely magnetic impurities even in the absence of Andreev coherence there is a 
peak on each side of zero energy. 

For Refs. there is no apparent electron-like component of the resonance in the DOS, 
thus local Andreev coherence is absent. In measurements of the LDOS near metal islands on 
BSCCO 01, however, both hole-like and electron-like peaks are apparent. This may indicate 
that the metal plays an important role in maintaining phase coherence at the surface, or that 
the metal overlayer is less disruptive to superconductivity than impurities in the plane. The 
presence of the in-plane Andreev processes, indicated by the electron-like peak, is essential 
to the operation of Josephson junctions. 

The on-site LDOS of Ref. is shown in Fig. |^(a). The second (hole-like) peak is not 
as clearly evident in the results of Refs. and thus may be peculiar either to the Zn 

impurity or to the impurity site in the BSCCO unit cell. Additional resonances around 
impurities can originate from additional orbital states around spatially-extended potentials 
or from spin-splitting near magnetic potentials. Note that an effective magnetic potential 
could also originate from a nonmagnetic impurity potential placed in a spin-polarized host 
electronic structure. 

Figure |^(abc) shows the best fit of dl/dV to the data of Ref. for phase incoherent 
Andreev processes and i a single-site nonmagnetic potential (1.375 eV), ii a nonmagnetic 
potential with onsite (0.360 eV) and nearest-neighbor (0.150 eV) values, and in a mixed 
nonmagnetic and magnetic potential (Vq = 0.825 eV, Vs = 0.550 eV). Also shown is the 
best fit using iv a pseudogap with no superconducting component and a mixed potential 
(Vo = 0.543 eV, Vs = 0.290 eV). The three panels show dl/dV (a) at the impurity site, (b) at 
the nearest-neighbor site along the gap nodes, and (c) along the gap maxima. Measurements 
of Ref. for (a) and (b) are shown; (c) is not available. 

The large size of the resonance on-site and the simultaneous presence of the gap feature 
occur because of junction normalization (equal resistance at —200 mV) and the finite width 
of the |0CTp (modeled as Gaussians of range O.SA, 3.8A, O.SA, and I.OA for i-iv). The very 
small LDOS in this energy range at the impurity site causes the tip to approach closer to 
the surface and (1) enlarge the apparent size of the resonance on-site, and (2) pick up the 
gap features from the nearest-neighbor sites. The relative size of the on-site resonance to the 
gap features is largely determined by the overlap of the nearest-neighbor Wannier functions 
with the tip when the tip is over the impurity site. 

Judging from the comparison with experiment, in and iv appear most in agreement, i 
does not have a second resonance, and whereas ii does show one in the proper location, its 
relative magnitude is incorrect. The smaller amplitude of the second resonance is obtained 
for Hi and iv because the overlaps with the STM tip are spin-dependent (|0|(x)p/|0|(x)p ~ 
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40). If the second peak were absent, either the impurity would lack magnetic character, or it 
would occur in less magnetic regions of the BSCCO unit cell. The remaining disagreement 
is in the amplitude of the resonance in (B), where iv is best, but still too small. 

Figure § shows the amplitude of the resonance as a function of distance from the impurity 
along the gap maxima (a) and the gap nodes (b). The squares are the data from Ref. 0, 
whereas the solid line corresponds to Hi, the dotted line to ii, and the dot-dashed line to iv. 
The plot for i looks identical to that of Hi. The agreement of Hi and iv are quite good along 
the maxima direction. The absence of a well-defined maximum at the nearest-neighbor in 
(a), which was pointed out in Ref. is due to the normalization procedure. An inset 
in Fig. ^ shows the difference between the junction normalized (solid) and unnormalized 
(dashed) dl/dV. The main discrepancy is with the amplitude of the signal along the node 
directions (b). 

Figure ^ shows the dl/dV of the resonance for iv; Hi is similar. The differences between 
the Figs. 2-4 and the measurements of Ref. |^ may be due to errors in the homogeneous 
electronic structure of BSCCO used in the calculation, particularly the low-energy electronic 
structure which dominates the longer-range LDOS. These errors may be due to inaccuracies 
in the model for the electronic structure measured by photoemission [|^, or they may be 
due to the neglect of other collective effects on the surface. Another likely source of error is 
that the electronic structure model of the host is not spin dependent. 

One of the possible mechanisms of a pseudogap is local AF order, such as occurs in a 
stripe. 1^ The magnetic component apparent in the impurity potential suggests this origin 
as well. If two nearby local moments are aligned parallel, then the resonances associated 
with them will hybridize and split |]T5|], whereas if they are antiparallel the resonances will 
be degenerate. A careful examination of the dl/dV for two Zn atoms near each other on 
the surface may clarify whether there is local AF order. 

The LDOS reported in Refs. [|]-|^ are best explained by the presence of a pseudogap state 
on the surface of BSCCO. The relative height of the electron-like and hole-like resonances in 
the DOS depends directly on the amplitude of local Andreev processes, and thus shows the 
degree of local superconducting coherence. This is of great practical interest, for the presence 
of these processes is essential to forming a proper Josephson junction across an interface. The 
information obtained about the superconducting state at the surface of BSCCO indicates 
the clear promise of future STM measurements near defects in other correlated electron 
systems. 

I would like to thank J. C. Davis for discussions and for providing the data of Ref. [Q. 
This work has been supported in part by ONR through contract No. N00014-99-1-0313. 
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FIGURES 




FIG. 1. Resonance energies for single-site (a) nonmagnetic and (b) magnetic impurity poten- 
tials. DOS for (c) nonmagnetic and (d) magnetic impurity potentials which produce a resonance 
at -1.5mV. Solid, 40 meV superconducting gap, dashed, 25 meV superconducting, 40 meV total 
gap, dot-dashed, 40 meV non-superconducting gap. (e) and (f), same as (c) and (d) except the 
dashed line corresponds to partial and the dot-dashed line to no phase coherence in the 40 meV 
superconducting system. 
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FIG. 2. dl/dV (a) on-site, (b) at the nearest-neighbor along the gap nodes, and (c) along the 
gap maxima, for impurity and host models i, ii, in, and iv. Data of Ref. is also shown. The 
order of vertical offsets (introduced for clarity) of the curves in (b) and (c) is the same as in (a). 
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FIG. 3. dl/dV as a function of position along (a) the gap maxima and (b) the gap nodes, ii 
is dotted, iii is solid, and iv is dot-dashed. Shown in inset is the difference between dl/dV when 
junction normalization is (solid) and is not (dashed) taken into account. 
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FIG. 4. Spatial structure of the dl/dV at the impurity resonance voltage (—1.5 mV) for iv. 
Horizontal is parallel to the gap maxima directions. 
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